Hole motion in the Ising antiferromagnet: an application of the recursion method 



Oleg A. StaryklH and George F. Reiter 
Texas Center for Superconductivity and Physics Department, University of Houston, Houston, TX 77204-5932 

(11 July 1995) 

We study hole motion in the Ising antiferromagnet using the recursion method. Using the 
retraceable path approximation we find the hole's Green's function as well as its wavefunction for 
arbitrary values of t/J z . The effect of small transverse interaction also is taken into account. Our 
results provide some additional insight into the self-consistent Born approximation. 
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I. INTRODUCTION 



The problem of the hole motion in a quantum antiferromagnet has become one of the central issue in developing 
the theory of high-T c superconductivity. The problem is an old one , but only recently has a good understanding 
of it been reached (for a review see, e.g., ||). The self-consistent Born approximation j|,U turns out to be extremely 
successful inpredicting the energy of the quasiparticle I^Q, mainly due to the vanishing of the low order vertex 
corrections |6|. Results obtained within this approximation agree well with exact diagonalization studies on small 
clusters B . The basic feature is that hole motion is strongly renormalized by the cloud of spin excitations (distortions) 
it causes, and results in a narrow (of the order of the superexchange constant J) band with minima at the (±7r/2, ±7r/2) 
points on the boundary of the magnetic Brilloin zone. More recent work by one of us Q has obtained the wave function 
of the hole within the same approximations and is also in agreement with exact diagonalization results [jTo[ . 

If the radius of the magnetic polaron is small enough, which happens at J ~ t, its wavefunction can be constructed 
variationally in direct space |pd| , |l2| , and the bandwidth becames of order t 2 / J . 

It was found long ago that hole motion in antiferromagnet with strong Ising anisotropy produces a string of the 
overturned spins, which confines hole to the origin of its path |||, making the hole's motion completely incoherent. 
The possibility of curing the spin background by the hole going one and a half times around an elementary loop on 
the lattice |13| does not change this picture significantly, as the effective mass of such coherent motion is extremely 
high. Surprisingly, the string picture seems to survive the limit of the isotropic Heisenberg model, where quantum 
• fluctuations that can cure overturned spins are most effective. This shows up in the subleading peaks above the 
£^ \ quasiparticle peak in the hole spectral function 

The subject of the present work is to study the hole motion in an antiferromagnet with strong Ising anisotropy 
using the recursion method. The method, also known as Lanczos technique, was developed initially for the electronic 
structure calculations of disordered systems |jl4|,[l5), and later generalized for the finite-lattice calculations of the 
strongly correlated systems HJl(|. We show that the method is particularly well suited for the problem at hand, and 
derive many of the known as well as some new results within its framework. To keep things as simple as possible 
we treat the t — J model within the linear spin wave approximation, which has no formal justification in the Ising 
limit. Nevertheless, we show that physics of the problem remains essentially unchanged by this drastic approximation. 
Moreover, the spin wave formalism permits us an important direct comparison with the results of the self-consistent 
Born approximation approach Q . 

The recursion method consists in the following juj. Given the Hamiltonian and the initial vector |1 >, it generates 
a new basis of vectors according to the rule 

\n+ 1 >= H\n > -a„\n > -b 2 r \n - 1 >, (1) 

with |0 >= 0, and b\ =< 1|1 >. The coefficients in the recurrence are calculated from 

_ < n\H\n > , 2 _ < n\n > 
< n\n > <n-l\n-l> 

Remarkably, in the new basis the Hamiltonian has a tridiagonal form 

H nn = a n , H n -x, n = i?„, n -i = b n , H nm = for all |n - m| > 1, (3) 
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and thus describes a fictitious semi- line of "atoms" with local orbitals |n >, diagonal energies a n , and hopping elements 
b n to the n+ 1 "atom" . Then the diagonal Green's function G u (w) =< — > takes the form of a continued 

fraction 

Gu(u) = hl h2 (4) 



u> — a-i ' 

We derive analytical expressions for the hole's Green's function and wavefunction in the next Section. The effect of 



small transverse interaction Jj_ is taken into account in Section III. We find that hole motion becomes coherent, and 



describe it for different limits of t/J z ratio. Comparison with previous works is done in the Conclusion. 



II. THE ISING LIMIT 



Let us consider the motion of the hole coupled to the localized spins on the lattice . In the linear spin- wave 
approximation the Hamiltonian of the t — J model reads (||],|7j 



H = zt^^ 4>(k, q) (c£c k - q a q + c+_ q c k a+^ + J z ^ w(fc)a^a fc , (5) 

k,q k 

(j>(k,q) = u q -y k - q +v q ^k, 

uj(k) = yjl- a 2 7 2, lq = -(cos(g x ) + cos{q y )) 

where we absorbed a factor zS in our definition of J {z is a number of nearest neighbours), u q and v q are the usual 
Bogolubov coefficients 



U " = V^T' v « = - S9n ^^)-> (6) 

and a = J±/J z . Each hopping of the hole produces emission or absorption of the spin excitations. 

For the problem of single hole motion in antiferromagnet it is natural to choose as a starting vector the state with 
one hole and no spin deviations, i.e. |1 >= c^\vac >, where \vac > denotes vacuum for both hole and magnon 
operators. Thus \vac > is simply the Neel state. The next vector is then 

|2 >= H\l >= H t \l >= ztY,4>M4- q a+\vac > (7) 

9 

Clearly, < 1\2 >= 0, and b\ = |fjf| = {ztf ~£ q <t?{p, «)• Acting by H on |2 > one finds 

H t \2 >= (zt) 2 ^2(t) 2 (p,q)c+\vac > +{ztf ^ <f>(p, q)<f>(p - q, qi)c+_ q _ qi a+ a+\vac > 
9 

Hj\2 >— ztJ u{q)4>{p, q)c+_ q a+\vac > (8) 

9 

a =j EMiW(p,q) 

Following the recipe described in the Introduction one obtains 

|3 >= {zt) 2 ^4>{p,q)4>{p- q,q{)c+_ q _ qi a^a+\vac > 



t T s^f J2 q UJ (q)4> 2 (p,q)\ ,, . , ,, 
+ztj^ iw(g) — ^ I Hp,qK- q 4\ vac> 



(9) 



The procedure becomes untractable very fast because of the branching: each action of H on a state with n-magnons 
produces again a state with n-magnons (Hj term) as well as states with (n+1)- and (n— l)-magnons (H t term), whose 
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coefficients have a complicated momentum dependence. The exception occurs when ui(q) is momentum-independent, 
i.e. in the Ising (a = 0) case. In that case one has 02 = J z and 




|3 >= (zt) 2 ^2(f)(p,q)(f)(p - <7,gi)c+_ 9 _ 9l a+a+|i;ac >, (10) 
9,91 

where now <f)(p, q) — Jp—q, and u q = 1, v q — 0. The norm of this state is 

< 3)3 >= {zt) A ^2 ^2 0(Pj 9)0(P - 9, gi)0(p, O0(P _ h h) < vac^ia^a^a^vac > 
9,91 Mi 

= {ztfY J ( t> 2 ip,d)<p' 1 {p-<i^) (11) 
9,91 

due to the property 52 7p-9i7p-<j- 9 i = 0- Thus, among two possible pairings of a-operators in jll| ) only the 
"diagonal" one (h = Qi, I = q) contributes. Note also that ~^2 qi lp- q - qi = §j and hence 63 = zt 2 . Continuing one 
finds 

= S,, gl Kg) + ^(gQ) 2 (p. g)^ 2 (p - g. gi) = 2J 
3 2 E,, gi 4?{p,qW{p- g,gi) 

|4 >= (zi) 3 ^ <^0-g-gi,g2)^0-g,gi)<?!'(p,g)Cp_ 9 _ 9l _ 92 a+a 4 " 1 a+|t;ac> (12) 

9,91,92 

Now, however, in addition to the "diagonal" pairing {U = q^) a new one (/ = q%, h = qi, h = 9) appears, for example, 

< 4|4 > = (zi) b <j (^) a + 7p-97p-92 7p-9l-92 7p-9-9l7p-g- gi -g 2 | (13) 

This pairing describes hole jumping around elementary loop (plaquette on a square lattice) 1.5 times jl3). The 
triple integral in ( |l3| ) is calculated to give (|) 3 (2 + cos^ — p y )) (2 + cos(p x +p y )). It describes an effective hole 
hopping along the same sublattice. The hopping is maximal for p — (0,0), and minimal for p = 0) or (0,7r). 
Propagation via closed loops was studied in great detail in Ref. Jl7j . For J z ^> t the effective hopping corresponds to 
the tunneling through a potential barrier, (it scales as exp [— 4.6( J^/i) 1 / 2 ] p7| ), and leads to the band with minima 
at p = (0,0). In the opposite limit (t ^S> J z ) hole can travel very far before noticing the confining potential due to 
the finite J z , and the weight of the closed- loop paths among the long self-retraceable paths is very small. This point 
of view is supported by the results of |lfi|] , where it was found that the ground state energy of the one-hole system 

scales as J z ^ 3 in a wide range of J z /t (1 > J z /t > 5 ■ 10~ 3 ), which is a characteristic feature of the retraceable-path 
approximation (see our discussion below). 

Based on these arguments we will omit all closed-loop contribution in the following. This corresponds to solving 
the problem on a Bethe lattice, or in the limit of infinite dimensionality fl9| . Diagrammatically this approximation 
consists in neglecting the vertex corrections Hf7J. 

With this simplification in mind the procedure can be iterated infinitely and results in the following expression for 
the coefficients of the recursion 

o-n — in — l)^z, b n — zt 2 , for any n > 2. (14) 

As should be clear by now, linear n-dependence of a n comes from the simple fact that Hj\n >— (n — 1) J z \n >. On 
the other hand the n-independence of b n is due to the particular form of the coupling j p — q . The strong momentum 
dependence of this coupling favors only one particular pairing out of (n— 1)! possible at the n-th step of the recursion 
procedure. This particular pairing is nothing but the retraceable path approximation of Brinkman and Rice [0: 
first the last excited magnon is absorbed, then the one before last, etc... If the hole-boson coupling is momentum 
independent, as happens, for example, in the Holstein model, the outlined procedure leads to b 2 n ~ n p0|) , and, as a 
result, to a different physics. 

So far we have found that Hamiltonian H has a simple tridiagonal form in the space spanned by the vectors \n > 
(see Eq.(^)). Then for any vector |\f >= ^n\n > of this space the Schrodinger equation H\^> >= > holds, 
or 

(E - a„)*„ - - 6„+i*„+i = 0, (15) 
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with boundary condition — (remember that |0 >= by construction). In small- J (large-n) limit we may use a 
continuum approximation to find 

±Vi*?rir±V2**» + (n-i) J z y„ = EV n (16) 

an 1 

Which is just the Schrodinger equation for a particle on a semi-line in the linear potential. Interestingly, index n plays 
a role of coordinate along the path of the hole - this correspondence is exact because we neglected all closed loops. 
The +(— ) sign in ( |l6| ) corresponds to choosing energy E close to the upper (lower) edge of the band. At the bottom 

2 /3 

of the band one finds J z -behaviour: 

E n = -Isfzt -J z + /? n jf/3(yit)V3 (17) 

here — n are the zeros of Airy function. It is known that the spin-wave approximation overestimates the bandwidth 
and more accurate treatment of the constraint leads to the replacement of \fz by the y/ z — 1 in the formulas 
above. 

The result obtained is by no means new and goes back to the sixties M. 

We would like to demonstrate a different approach here. Consider the one-hole Green's function, which has a 
continued fraction form (see Eq.(|4|)). 

Gn («) = X — (18) 

u ; TP — 

A little thinking shows that it is equivalent to 

G{uj){uo- zt 2 G{u- J,)) = 1 (19) 

which coincides with the self-consistent Born approximation ||, in agreement with our expectations. Let us try the 
ansatz [^lj 

y/Zt T(CJ + J z ) 

Then Eq.([l9|) takes the form of the difference equation 

T(w- J z ) + T{uj + J t ) = — ^-T( W ), (21) 

\/zt 

which is a well-known recursion relation for the Bessel functions [p2[ , and hence T(w) = AJ^^l (j^j^j- This 
immediately gives for the Green's function 



J. 



( 2y7t 

ft V J. 



G u (u) = Y' \ (22) 

This form is a direct consequence of the linear dependence of a n on n (see Ref. ]23|), which, in turn, is an intrinsic 
feature of the retraceable path approximation. A similar expression was obtained within the retraceable-path approx- 
imation (without using spin- wave transformation) in Ref. p4j . Notice that energy appears only in the indexes of the 
Bessel functions. 

Let us look on the -j- ^> 1 limit first. Here we need the double asymptotic expansion of the Bessel function at large 
value of argument and index [B2l 



./^ cos <?(/«) ~ y ^ sin q cos {^(sinfl - 9 cos 9) - ^ j (23) 
H > 0, < 6 < tt/2 



Somewhat lengthly calculations give an expected answer 
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oo 

G "^ = Wt^—, (24) 

v i/=0 

uj v = -2y/zt ~J Z + p v Jf /3 (Vi*) 1/3 , 

in complete agreement with results of the Schrodinger equation (p"7j). Notice that residue of the poles is ^-independent, 
and is given by ^j^, supporting the "dominant-pole" approximation of We note in passing that the quasiclassical 

expression for the zeros of the Airy function j3 u = (^-{y + |)) 2 ^ 3 works extremely well even for v = 0: it gives ~ 2.32 
compared to the exact value /3q = 2.34. 

In the opposite -f < 1 limit one can use standard small argument expansion 

to find the perturbation theory result Gn(w) = (u + 

Recently one of us || derived a wavefunction that corresponds exactly to the self-consistent Born approximation, 

>= Zk (c+ + N- 1 ' 2 l)Gk- q {e k - co q )c+_ q a+ + ... (26) 

+N -n/2 (zt) n ct>(k,q)G k - q (e k -u} g )...<f>{k-q-qi- ... -q n - 2) q n -i) 

q,qi, ...q n -l 

xGt-„i-...-fc-i( f t -Q) q - ... - u qn _Jc+_ q _ qi _ _ qn _ iC L+ ...a+^Jlvac > . 

Here e k = Sfe(efc) is the quasiparticle energy, and is the self-energy in the SCBA. The quasiparticle spectral weight 
is given by (zfe) 2 = ^1 — ^%£^U=«k) ■ Clearly, in our case the role of the momentum k is played by the index v. 
The specific form of the Green's function ( p2f ) immediately gives us that the chain product of the Green's function in 
the last equation has a simple form G(u - J z )G{lu - 2J z )...G(u - n.J z ) = (-^)"J J j h£l)/j_ J ^(h£L). Using 
again asymptotic expansion of the Bessel functions we find for the hole in lowest state v = (t/J z 1) 



'I cba >~ - ^ 1/2 E^7,- 9 4_ q aj + 

_ + N -n/2 J- (-l)"W 2 n 7fe _ g ... 

q,qi,...q n -i 

x 7fc-,-...- 9 ,-i 4_ q _ qi ___ qn _ i a+...a+ n _ 1 )\vac>, (27) 



for n < (^t/J z ) 1/3 . 



III. SMALL J± LIMIT 

A virtue of the present approach is that we can include the effect of small transverse interaction J± , J± = aJ z <C J z ■ 
As was already discussed by analytical consideration is possible if one restricts to linear in a accuracy. Then 
ui(q) = 1 + 0(a 2 ), u q = 1 + 0(a 2 ), but v q = —\u^ q . Hence 

4>(P, l) = lp-q - + °( a2 ) ( 28 ) 

One finds for first a- and 6-coefficicnts 

ai =0, a 2 = J z , b\ = (zt) 2 ^ (ip-q ~ a7 P 7g7p-g) = zt 2 (\ - aj 2 ) (29) 

9 

The key property 7p-g7fc_„ = insures that the retraceable path approximation is exact at low order, i.e. that 
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J2 <t>{p,<a)<t>(p-ii, q)<t>(p-q, <n)<t>(p, q) = o + o( a 2 ). (30) 

Moreover, one can find 

/„ = ^2 & 2 (Piq)4> 2 (v~ q,qi)--4 2 (p- q- q\ - - qn-2,q n -i) 

g,gi,...g„_i 

= X! ( l ,2 (p,q)4' 2 (p- q,qi)--<t> 2 (p- q- qi - qn-^qn-2) 

q,qi,...q„- 2 

x X! ~ q ~ qi ~ ■■■ ~ q,i - 2 > qn -^ 



n-2 



z z \ z I ' H 



As a result 



- azJ5) + O(a^), n > 3 (32) 

A _ 

" 64 
q 

At the same time there are no changes in a n , a n = (n—l)J z . Due to the difference between b\ and &^ 2 ' ^ ne continued 
fraction expression for the Green's function is now equivalent to the following system of two equations: 

Gll{Uj) = a;-z*2(l-Q 7 2)G 22 M 

GmM = ^— (33) 

We introduced P — t 2 (l — azB) here. Comparing the last equation of j3^ ) with Eqs.( |l9| , p2"|) of the previous Section 
one finds the answer 

J U ,_ J (hHl\ 

G 22 H = --L ZIZEZy±J (34 ) 

In the J 2 <C i limit expression ([54j) takes familiar form 

= ^=£; —^r- (35) 

We want to consider the lowest hole state, and thus may approximate 

^ W - LJ V UJ-UJ ^ LUq-UJv 

The last sum can be approximated by the integral from to some v m , which, in turn, is determined by the condition 
Co Um ~ 1\[zt (upper edge of the incoherent band at J z = 0). In this way we find that last sum is given by — We 
substitute the result just found for Giiiio) into the first of the equations (|33]), and find after some algebra 

1 



4(tt - 2) y/zt uj-Cj 

Co = -2V^+ 2.34 J 2 J\V^) 1/3 ~ 2(l J ~2) (1 ~ ^ (3?) 



G 



The pole at tu describes a narrow, coherent band of bandwidth ~ 1.38J_l with minima at j p — 0, i.e. along the 
boundary of the magnetic Brilloin zone ||. Note that residue of the pole remains unaffected. It is a matter of short 
calculation to find out that in the opposite limit (J z ^> t) the answer is 

GnH= -rfs4 iT (38) 

The bandwidth is much smaller, zt 2 Jj_/ J 2 , and goes to the correct perturbative answer ~ zt 2 / J at the isotropic point 
Jj_ = J z (a = 1). 



IV. CONCLUSIONS 



A simple recursion method was used to calculate a single particle Green's function for the hole moving in the 
antiferromagnet with strong Ising anisotropy. By neglecting closed-loop contribution we reduced the problem to an 
exactly solvable one. In the Ising limit we find, in agreement with known results, that hole is confined to the origin of 
its path by the effective linear potential due to the overturned spins. We then take into account the small transverse 
interaction between spins and find that the hole motion becames coherent. We calculate the hole's spectrum and 
Green's function to the first order in Jj_/J 2 . This calculation explains the "dominant-pole" approximation of ||. We 
argue that the spin- wave approximation, employed in this paper, does not affect the essence of the problem, as can be 
seen, for example, from the comparison of our expression (23) for the hole's Green's function with formula (4) of Ref. 

f. Overall, our approach can be considered as an another way to arrive at the self-consistent Born approximation 
at least in the case of strong Ising anisotropy. Our results become exact in the limit of infinite dimension. 
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A.L.Chernyshev for valuable discussions and comments, and to V.Brazhnikov for help. 
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